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Abstract 

O 
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We present a simple method for the construction of exact ground states of 
generalized Hubbard models in arbitrary dimensions. This method is used 
to derive rigorous criteria for the stability of various ground state types, like 
the 77-pairing state, or Neel and ferromagnetic states. Although the approach 
presented here is much simpler than the ones commonly used, it yields better 
O ■ bounds for the region of stability. 
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In recent years there has been an increasing interest in analyzing systems of correlated 
fermions by constructing the ground state explicitly in certain parameter regimes, but ar- 
bitrary dimensions |L|- pi. In general those ground states have a rather simple structure. 
Typically one can speak of quasi-classical states where quantum fluctuations are not im- 
portant and there are no finite-size corrections to the ground state energy. E.g. the fully 
polarized ferromagnetic state has been in the center of attention. In |7j the corresponding 
program for superconducting states of the //-pairing type 0,0] has been initiated. These 
^-pairing states show Off-diagonal-Long-Range-Order (ODLRO) [[K]] which in turn implies 
the Meissner effect and flux quantization P~C|— |T2|| , i.e. superconductivity. 

The most popular method for the construction of exact ground states was introduced 
by Brandt and Giesekus |] and generalized by Strack and Vollhardt |§. The basic idea is 
to start with a Hamiltonian Tio with a known ground state \iPn) and add operators V^V 
where V annihilates | iPn)- Thus 7i + V^V also has | ip^) as ground state. It is obvious 
that for a given Hamiltonian 7i and ground state | ipN ) it is usually not easy to find an 
appropriate 7i and V . Application of this method requires a lot of guesswork and is not 
very systematic. In addition, one does not gain much insight into the underlying physics by 
this sort of construction. 

Ovchinnikov [^],^[ has used a different method to improve some of the results obtained 



previously. His approach is based on Gerschgorin's Theorem (see e.g. [13]) which gives a 
lower bound on the ground state energy of the Hamiltonian and thus complements the usual 
variational principle which gives upper bounds. If those bounds coincide each eigenstate 
with the corresponding energy is a ground state. 

In this paper we apply a much simpler and clearer method for the construction of exact 
ground states which we call Optimal Ground State Approach (OGS approach). The basic 
idea is to diagonalize the local interaction and make all local eigenstates which are needed for 
the construction of the global ground state | local ground states. This usually implies 
some restrictions on the interaction parameters, typically in the form of inequalities. We 
will show on a few explicit examples that this method is not only much simpler than the 
ones described above, but also reproduces or improves the results for all cases considered 
previously. We like to point out that the OGS method does not only work efficiently for 
states | ipN ) which are simple tensor products, but also for more complicated ones, like the 



77-pairing states or so-called matrix-product ground states of certain spin-1 chains JT^]. A 
more detailed exposition of the OGS approach and a thorough comparision with the other 
methods will be presented in a future publication . 

For any Hamiltonian 7i = hji 011 a lattice with L sites ((jl) denotes neighboring 
sites) it is usually quite easy to diagonalize the local interaction hj\. Suppose now that the 
lowest eigenvalue of hji is (e.g. by adding a suitable constant). Clearly is a lower bound 
for the global ground state energy. The global ground state [ Vjv ) is called optimal |L6| iff 

^Pn ) = 0, i.e. the ground state energy is just the lower bound found by diagonalizing the 
local interaction. 

In the following we will be interested in a generalized Hubbard model with 
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+V(n j - l)(n, - 1) + y (cj^q^ + cj^c^c^) 

+ J -f { s l s i + s f s l) + J'Sjs; + 1 K- + ni ) . (1) 

Here Cj CT and cj CT are the canonical Fermi operators, rij a = Cj a Cj a and rij = + riji are the 
corresponding number operators and the SU (2) spin operators are given by Sj = riji), 

S I = c li c iT and s t = c h c n- 

The first term in (fll) is the single-particle hopping, the second one is known as bond- 
charge interaction. U and V denote the on-site and nearest-neighbor Coulomb interaction, 
respectively. In addition, we included a XXZ-type spin interaction with exchange constants 
J xy and J z between nearest neighbour sites and a pair-hopping term Y . \i is the chemical 
potential and Z the coordination number of the <i-dimensional lattice. 

The local Hamiltonian can easily be diagonalized. Denoting an empty site by 0, a site 
occupied by an electron with spin a ="[*, j by a, and a doubly occupied site by 2 we find the 
following 16 eigenstates and their respective energies: 

state eigenvalue 
100) E 1 = IL + v 

\a0) ± |0a) 4 ±) ==F* + | 

I \ w U J 2fi 

ITD + Ilt) ^ = "|? + if "T + Y 

|20)-|02) E 5 = ^-V-Y+^ 
|V±> = « ± (|U)-||t)) + (|20) + |02)) 4 ±} 

|a2)±|2a) ^ = ±(t - 2X) + % 

122) E 8 = — + V + ^ 
17 8 2Z Z 



with 



and 



(2) 
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— — ±V — Y — — ^ ±23 

a± = z - 2 - 4 P (4) 

* 4(X - t) { V 



where /? = ^ - g - f - ^ - i) 2 + 4(X - t)2. 

For the special case t — X, which will of some importance in the following, this simplifies 

to 

hM = IU>-IIT>, 

|^) = |20) + |02), (5) 

with corresponding energies 

„(+) _ _U_ Jxy Jz 2/x 
^ ~ 2Z 2 4 + Z ' 

& = *z- v + Y + %- < 6 » 

First of all we look at the ^-pairing states with momentum P, 



l^> = WlO), Vt-E^A- (7) 

3=1 

As already mentioned these states exhibit ODLRO and are thus superconducting [[J. 

The ^-pairing state with P = is an eigenstate of 7i for £ = X, 2V — Y [j/j. In order 
to make it an optimal ground state we first observe that IV'jv ) can be built completely from 
the local states 1 00) , 1 22) and 1 20) + 1 02) . All three states are already local eigenstates with 
energy E = E\ = E$ — E s — U/2Z ± V, if we choose fi = 0. Demanding that E is 
the local ground state energy one recovers the result obtained in |7| by using Gerschgorin's 
theorem: 

V < 0, 

-| > max(2|t| + 2V, V - |, V + ^ ± (8) 

The state with momentum it is an eig enstate of H for Y = -2V. is now built from 

the local states 1 00) , 1 22) and 1 20) — 1 02) . These states have to be made local ground states. 
Again for \i = they all already have the same energy E = E x = E 5 = E 8 = U/2Z ± V. 
All other energies have to be larger. This yields the following inequalities: 

V < 0, 

-- > max(2|t| + 2V, V - ^, 2\t - 2X\ ± 2V, V - % + ^, 

7, ~ V 1 1 4 ' 1 1 9. A 
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J*y Jz_ it-*)* } (g) 

2 4 V 1 
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There exist also ^-pairing states with momentum P ^ 0, ir. These appear for instance 
as ground states of model ([l]) where t = X and U < — 4t, but all other interaction constants 
are zero |7|, ^7| , P^] . These states are eigenstates only for t = X and Y = V = 0. It is 
straightforward to derive a similar inequality as for the states with momentum or it. It is 
just (|J) with V put equal to 0. 

There are other interesting states for which the OGS method can be applied, e.g. the 
following paramagnetic, Neel, charge-density-wave and ferromagnetic states at half-filling 
(N = L, where N = Ylf=i n j is the total number of particles), 

\para) = IJ 4t II c};|0), (10) 
jeA jeA' 

|iVeeO = IIs t TlIS t |l )' ( n ) 

jeB jeB' 

\CDW) = l[c] ] c] l \0), (12) 

jeB 



= n^tio)- (13) 



A and A' are arbitrary disjoint sets of lattice points which together span the whole lattice. 
The states | Neel) and \CDW) are defined on a bipartite lattice with odd and even sublattices 
B&ndB'. 

If we want \para) to become an optimal ground state we have to make all the states \crcr') 
with a, a' =1,1 local ground states. From (§) we see that we have to choose J xy = J z = 
and t = X in order to get the correct local ground states. These then have the energy 
E 3 = E 4 = E { 6 +) = -U/2Z + 2(i/Z. All other energies have to be higher which finally leads 
to the condition 

Y\-V + $f). (14) 

For Y = and V > this problem has already been investigated in 0. In this case the 
bound (0) is better than the one found in 0. In fact, for /i = it is exactly the improved 
bound found by Ovchinnikov || using the Gerschgorin approach. 

Similarly, for the Neel state \Neel) we get the restrictions t = X and J xy = (from 
E A = Eq + ^) and the inequalities 





2|f|-£ + %U + |y|-£). (is) 

Turning to the state \CDW) we see that this state can be constructed completely from 
the local states 1 20) and |02). In order to make these local ground states we have to choose 
t = X, but also Y = 0. The local ground state energy is then E 5 = E^ = U/Z — V + 2/j/Z. 
The condition that the remaining energies are higher leads to the inequalities 
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(16) 

which reduces to Ovchinnikov's improved result |J for J xy = J z = 0. 

The OGS approach shows why a generalization to the non-half-filled case is difficult. 
Away from half-filling also other states (like |er0) ± |0c)) would have to be used in the 
construction of the global ground state. Thus the local ground state degeneracy would have 
to be larger. In most cases this leads to too much restrictions on the interaction parameters. 

For the fully polarized ferromagnetic state (at half-filling) \F) the local ground state is 
only \o~o~). Therefore we do not get restrictions like J xy = J z = or t = X in this case. The 
local ground state energy is E 3 = —U/2Z + J z /4 + 2fi/Z. This is a lower bound for the 
other local energies, leading to the inequalities 

Jz ^ I Jxy I > 

-> m ^2\t\ + - + Y ,2\t-2X\ + -- Y ,- [ -V + —,V + Y + -, 



If one is interested in the sector with a fixed particle number N, one can regard the bounds 
as function of /i and try to find the value of \x which optimizes these bounds. If a state is a 
ground state of 7i at some fixed particle number N, then it is also a ground state of 7i + fiN . 
Hence we only need to require that a state is a ground state for some value of fi, the result 
for arbitrary fi follows then immediately. Therefore one can sometimes improve the results 
of the optimal ground state approach, by first finding inequalities for an arbitrary value of 
the chemical potential, and subsequently optimize with respect to /i. For example, if we find 
inequalities U > a + \i and U > b — ji, then the best value of \i is (b — a)/2, and thus the 
inequality U > (a + b)/2. 

For (|14D , (|15D , and ([U]) one obviously has to choose [i = 0. In the case of the ferromag- 
netic state the second inequality of ( |T7| ) can be replaced by 

J z , nvl V 3J Z . . V 3J Z 

- 2X\ + — , t - 2X\ + — , t + — , 

1 2' 1 1 2 8' M 2 8 

T-^ + y + T^- y + T"TO)' (18) 

This bound improves the one found in Q]. 

We also applied the Gerschgorin approach to the states considered here. In all cases the 
results are not better than the ones obtained from the OGS approach with optimization 




with respect to fi fll5| . 



In summary, we have shown that the OGS approach is - dispite its simplicity - a very 
useful tool for the construction of exact ground states. It allows to reproduce or even im- 
prove bounds for stability regions previously found by different methods and yields better 
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insights into the physics. E.g. if the interaction parameters are changed from a situation 
where the inequalities are satisfied to one where they are not, the ground state typically first 
is non-degenerate or doubly degenerate (except \para), which is always a highly degenerate 
ground state), becomes highly degenerate at the OGS bounds (where levels of the Hamilto- 
nian cross), and finally cease to be the exact ground state. In addition, by examining which 
local eigenstates become the new local ground states one can try to predict some properties 
of the new ground state once the bounds are violated. We believe that these bounds are as 
sharp as possible. If an inequality is violated we have a new local ground state which might 
be used to construct a variational state with an energy lower than that of the "optimal" 
state. We will come back to these points in a future publication |TJ| where we also discuss 



why the OGS approach usually yields bounds which are as least as good as those of the 
Gerschgorin approach. 
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